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Abstract 

The Cosmic Infrared Background (CIB) provides an opportunity to constrain many properties of 
the high rcdshift (z > 6) stellar population as a whole. This background, specifically from 1 to 200 
microns, should contain information about the era of rcionization and the stars that are responsible for 
these ionizing photons. In this paper, we look at the fractional anisotropy (SI /I) of this high redshift 
population, where SI is the ratio of the magnitude of the fluctuations and / is the mean intensity. We 
show that this can be used to constrain the escape fraction of the population as a whole, because the 
magnitude of the fluctuations of the CIB depends on the escape fraction, while the mean intensity 
does not. This results in lower values of the escape fraction producing higher values of the fractional 
anisotropy. This difference is predicted to be larger at longer wavelengths bands (above 10 microns), 
albeit it is also much harder to observe in that range. We show that the fractional anisotropy can also 
be used to separate a dusty from a dust-free population. Finally, we discuss the constraints provided 
by current observations on the CIB fractional anisotropy. 

Subject headings: cosmology: theory — diffuse radiation — galaxies: high-redshift — infrared: galaxies 



1. INTRODUCTION 

Modern cosmology is now able to constrain details 
about the era of reionization. Observations show that 
the reionization of the universe occurred early and was 
extended in time, with an equivalent of an instant aneous 
reionization at z ~ 11 l)Komatsu et al . 2009, 20111). Stars 
arc a likely candidate for being responsible for the major- 
ity of rcionization because they are efficient producers of 
ultraviolet photons. Thanks to a wave of modern, sensi- 
tive telescopes, we can begin to observe and understand 
the frontier of rcionization, along with these high rcdshift 
stellar populations (z > 6). For example, we can observe 
these galaxies directly via high redshift surveys, which 
can now routinely identify a population of bright galaxies 
up to a redshift of about z ~ 8. However, these surveys 
can only locate those galaxies that are both above the 
limiting magnitude and common enoug h to be present 
in th e survey field. It is now thought (iBouwens et al.l 
l20Ul IRobertson et~aTl f20Tol : iFernandez fc Shulll 120 111) 
that reionization needed a large population of smaller 
galaxies below the current detection limits. Because we 
cannot yet observe these galaxies directly, we can in- 
stead look for their cumulative light, which should exist 
as background radiation. Because reionization is said to 
have occurred around z ~ 11, the photons responsible 
for reionization should be present in the Cosmic Infrared 
Background (CIB). The spectral peak of this radiation 
is around the Lyman— a line, which will be redshifted to 
1 — 4 microns. However, continuum emission will create 
an extended tail at longer wavelengths. 

Here, we discuss the CIB from about 1 to 200 microns. 
The majority of the CIB will be emission from sources be- 
low z ~ 6, such as our Galaxy, foreground galaxies, and 
other sources of infrared light, such as zodiacal light. If 
these sources can be subtracted away to a high precision, 
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it is possible that the remainder could be from the era of 
rcionization, and if so, could tell us about the properties 
of these high redshift stars. 

There have been many attempts to theoreti- 
cally model the high rcdshift component of the 
CIB, especially in the near-infrared, fro m the 
dSantos et all 120021: iMagliocchetti et al 
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In this 



paper, we ex- 
amine another way to analyze the CIB, the fractional 
anisotropy, which is the ratio of the fluctuations to the 
mean. By looking at the fractional anisotropy, many 
free parameters arc removed and more information 
about this elusive stellar population can be extracted. 
Specifically, we discuss using the CIB as a probe for the 
escape fraction of ionizing photons. Finding the escape 
fraction is important for understanding reionization 
and its duration. There have been several attempts to 
measure the escape fraction through analytical models, 
simula tions, and observations (see IFernandez fc ShuTl 
(|2011h and references within). These papers have shown 
that the escape fraction appears to vary greatly from 
galaxy to galaxy. Therefore, instead of trying to measure 
the escape fraction of an individual galaxy, here we 
discuss the average escape fraction of all galaxies, which 
will give more of a global view of rcionization. In ad- 
dition, the fractional anisotropy can reveal information 
about the dust content of galaxies, which is mostly 
unknown at high redshifts. 

We describe our simulations in section [2] and our mod- 
els in section [3l In section 21 we discuss our method for 
finding the mean CIB, the fluctuations of the CIB, and 
the fractional anisotropy. In section [5j we discuss our 
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results of the fractional anisotropy for various bands. 
In section [6j we discuss the most recent observations. 
We conclude in section [7] Throughout this paper, we 
use the cosmological parameters (fi m , fib, /i)=(0.27, 
0.73, 0.044, 0.7), consistent with the simulations from 



llliev et all (pOia ). which are based on th e WMAP 5-year 
resul ts and other available constraints (|Komatsu et al.l 
l2009h . 

2. THE SIMULATIONS 

In order to predict the angul ar power spectrum of the 
CIB, we used simulations from llliev et al.l (|2012f ). which 
arc N-body simulations combined with radiative trans- 
fer, which allow us to sec how sources are affected by 
the reionization process. The high resolution of these 
simulations (with a minimum mass of lO 8 Af0) allow 
us to also include Jeans-mass filtering on low mass ha- 
los. This effectively allows suppression of star formation 
within small halos (10 s — 10 9 M & ) because of elevated 
gas temperatures that could be caused by the proximity 
to other star forming galaxies. These simulations have 
a box size of either 114h~ 1 Mpc (for cases with suppres- 
sion of small sources) or 37h~ 1 Mpc (with no suppression, 
where M min = 10 8 M Q , or complete suppression, where 
M m i n = 10 9 Mm). These simulations are summa rized in 
Table HI (llliev et all 12071 iFernandez et alJl2012j ). 

To describe the stellar populations within these halos 
and their relationship with their environment, we define 
a parameter / 7 , which describes the number of ionizing 
photons produced per stellar atom that can escape the 
galaxy and reionize the intergalactic medium (IGM). / 7 
is defined as a product of the star formation efficiency, or 
the fraction of baryons that are in stars (/*), the escape 
fraction (/ esc ), and the number of ionizing photons per 
stellar atom (Ni): 



f'y fescf*Ni. 



(1) 



We allow / 7 to have different values, dependent on the 
mass of the halo, assuring it is consistent with reioniza- 
tion. 

3. OUR MODELS 

True first generation stars are metal free (Population 
III stars). As time goes on, stars die and enrich the uni- 
verse, and eventually, these Population III stars give way 
to stars with metals (Population II stars). It is unclear 
when this happens, and this process is probably very in- 
homogencous. Therefore, we assume two limiting cases - 
all of the stars from 6 < z < 30 are either Population III 
(Z = 0) or Population II (Z = 1/5OZ ) stars. 

In addition to the mctallicity, there is uncertainty for 
the mass spectrum of these stars. These stars could be 
very large, or they could be similar in size to what we see 
today. To model these two e xtremes, we ch oose either a 
heavy, Larson mass function ([Larson! 1 1 9981 ): 



/(m) oc m 1 ( I + 



-1.35 



(2) 



with mass limits of mi = 3Mq, mi = 500Mq, and m c = 
250M Q to mod el a population of large stars, or a Salpctcr 
mass function (|Sarpeterlll955f ): 



with mass limits of mi = 3Mq and 1112 = 15OM0 to 
simulate a mass spectrum similar to what we see in the 
local universe. 

If we combine our limiting cases for both mass and 
mctallicity, we can establish our two limiting stellar mod- 
els: Population III stars with a Larson mass spectra, and 
Population II stars with a Salpeter mass spectra. In re- 
ality, these stellar limits are extreme. In addition, we 
expect stellar properties to be inhomogeneous through- 
out redshift. However, these examples were chosen as 
limiting cases: which represent a population with the 
smallest and largest amplitude for the angular power 
sp ectrum of a larg e rang e of models, studied in detail 
in IFernandez et aLl (|2010j ). We would expect the actual 
amplitude for the angular power spectrum to lie between 
these extremes. These populations are summarized in 
Tabled 

4. COMPUTING THE FRACTIONAL ANISOTROPY 
4.1. The Mean Cosmic Infrared Background 

Now that we have our stellar and galactic models, we 
are in a position to calculate the fractional anisotropy, 
SI/I. To do this, we must compute both the mean in- 
tensity and the angular power spectrum of the CIB. The 
mean CIB is a combination of emission from the star, 
which is modeled as a stellar blackbody, and emission 
from the nebula, which is a combination of the Lyman— a 
line, two-photon, free-free and free-bound emission, and 
with the possibility of emission from dust (see section 
14. 2[) . This nebular emission is either produced within 
the halo itself, or within the IGM if some fraction of the 
ionizing radiation (f eS c) escapes the halo. However, the 
mean CIB does not depend on / esc , since the nebular 
emission is the same, regardless of whether it is from the 
halo or the IGM. 

Each emission process (stellar, free-free, free-bound, 
two-photon, and the Lyman— a) was mode l ed an alyti- 
cally (for details, see IFernandez k, Komatsul (|2006f )) and 
integrated over a range of rcdshifts from 6 < z < 30. 
The total intensity (I) is then: 



I=±(f* 

4tt V n m 



-halo 

:Pm 



z) 



H( z )(i + Z y 

x [l*(z) + P f (z) + V\z) + P^z) + ^<*{z)] . (4) 

Here, p\f°(z) is the mean mass density collapsed into 
halos from the simulation. The luminosity per stellar 
mass, l(z), is given for each component of the luminosity, 
* for stellar, // for free-free, fb for free-bound, 27 for 
two-photon, and Lya for Lyman-alpha emission. The 
luminosity of any component " a" can be written asQ : 



TO = 



dlnp,(z) dmf(m)L^(m)T(m) 



dt 



Imi dmf{m)m 



(5) 



/(m) oc m 



-2.35 



(3) 



The luminosity of each emission component (£"), and 
the stellar lifetime (r(m)) are integrated over a mass 
spectrum of stars. The first part of this expression is 
the inverse of the star formation time scale, tsF( z ) = 

1 This expression is only valid if the average stellar lifetime is 
always less than the star formation time scale tgp, which is true 
for the cases we are c oncerned with. For more information, see 
IFernandez et al.l 120TOT ). 
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Simulation Name 


Box Size 
(Mpc) 


Minimum 

Halo Mass (Mq) 


Suppression 


J ~y , large 


J "i , small 




f 


Partial Suppression, 
















High Efficiency 


163 


10 8 


Yes 


10 


150 


8.3 


0.080 


Partial Suppression, 
















Low Efficiency 


163 


10 8 


Yes 


2 


10 


6.7 


0.058 


No Suppression 


53 


10 8 


No 


0.4 


6 


8.6 


0.078 


Complete Suppression 


53 


10 9 


Yes - complete 


12 





8.3 


0.071 



TABLE 1 

Radiative transfer simulations used in this work. z ov is the redshift of overlap, where reionization is complete, and t is 

THE ELECTRON SCATTERING OPTICAL DEPTH, f^. small Is FOR HALOS THAT ARE BETWEEN 10 8 AND 1O 9 M0 WHILE f^. large IS FOR HALOS 

ABOVE 1O 9 M . 





/esc 


/» - Pop II Salpeter 


/» - Pop III Larson 


10 


0.1 


3.8 X 10~ 2 


4.0 x 10~ 3 


10 


0.3 


1.3 x 10" 2 


1.3 x 10" 3 


10 


0.5 


7.7 x 10~ 3 


8.0 x 10" 4 


10 


1 


3.8 x 10~ 3 


4.0 x 10" 4 


150 


0.1 


5.8 x lO" 1 


6.0 x 10~ 2 


150 


0.3 


1.9 x lO" 1 


2.0 x 10~ 2 


150 


0.5 


1.2 x 10 _1 


1.2 x 10 -2 


150 


1 


5.8 x 10~ 2 


6.0 x 10~ 3 



TABLE 2 

THE PROPERTIES OF THE STELLAR POPULATIONS. /* WAS SET TO BE CONSISTENT WITH REIONIZATION. 



[ rfln P»( z ) ]-i _ This star formation time scale is unknown, 
but we assume a value o f 11.5 Myr , consi stent with the 
value from simulations of llliev et al.l (|2012D . This expres- 
sion then reduces to: 



i dm f( m ) L »( m ) T ( m ) 



(0) 



t SF(z) J mi 2 dmf(m)m 
(|Fernandez et al.H2010fl Pl. 

4.2. Dust 

We do not know how much dust exists in high red- 
shift galaxies. Molecular gas is already observed at 
z ~ 5, an indication t hat dust is present at those redshifts 
(jRiechers et al.ll2010f ). Because dust will affect the spec- 
tra of high redshift galaxies, the fractional anisotropy 
may also change. In order to see if our results are affected 
by dust, we compute the spectra expected if the radiation 
field is further reprocessed by a dusty medium. We gen- 
erate d a dust spectrum using DustEM (|Compiegne et all 
120111) predicted for a galaxy with a high metallicity 
and m inimal destruction of dust grains (jCompiegne et all 
2010). In reality, the low metallicity and hard radiation 
fields expected at high redshift will lead to a dust contri- 
bution that is less than the one modeled here. In addi- 
tion, DustEM models are computed in the optically thin 
limit, so therefore, the dusty SED we obtain is the upper 
limit for the amount that dust will redden. Our dusty 
model represents the extreme model for a dusty galaxy, 
with our case with no dust representing the opposing 
limit. 

2 The value chosen for tgp will change the amplitude of the 
luminosity, which will both affect the mean and the fluctuations of 
the CIB. A benefit of taking the fractional anisotropy is that the 
dependence on tgp will nearly cancel out. For more information 
on the dependence of the luminosity on tgp, see section 6.1 of 
IFernandez et al.l EoTOl') . 



4.3. Fluctuations in the Cosmic Infrared Background 

The next step is to compute the angular power spec- 
trum. These fluctuations will arise from both the emit- 
ting halos and their surrounding HII region s within the 
IGM. As shown in IFernandez et al.l (|2010tl . the fluctu- 
ations from the IGM are probably quite small (from 2 
to 7 orders of magnitude smaller than that of the halos 
themselves) so can safely be ignored. 

The angular power spectrum Ci can then be written 
as: 



Ci 



(4tt) 2 



$ 1 



dz 



H(z)r 2 (z){l + z) 4 



x [&°(z){l*(z) + (l-f esc )L(z)}Y 



xb 



eff 



r{z)'' 



Z Plin k 



r(z)' 



(7) 



Here 6 e // is the effective bias, Pn n is the linear matter 
power spectra, r(z) = c J* dz' /H(z') is the comoving 
distance, and the luminosity is: 

L(z) = l ff (z) + F b (z) + P\z) + ^(z). (8) 

The simulations provide both the halo bias and the linear 
matter density fluctuations. Note that the angular power 
spectrum depends on f esc . The angular power spectrum 
of these simulations was computed in Fernandez et al. 
2012. 

4.4. The Fractional Anisotropy 

The fractional anisotropy of the CIB is obtained by di- 
viding the angular power spectrum C; (shown in equation 
\7} by the mean intensity / (equation [J): 



51 1 1 = ^l(l + l)Ci/(2irP). 



(9) 
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Most of the free parameters then cancel out, including 
the star formation efficiency /». The luminosity l a , how- 
ever, will only cancel out when f esc = 0. Therefore, the 
fractional anisotropy serves as a test to constrain f esc . 

5. RESULTS 

At large values of I, the minimum mass of the star 
forming halos and suppression history will change the 
shape of th e angular power spectr um due to non-linear 
bias effects (jFernandez et al.ll2012l ). Since the minimum 
mass of these star forming halos is unknown, we com- 
pute the fractional anisotropy for I = 3000, avoiding any 
flattening or steepening of the angular power spectrum 
that could occur at larger 1. 

The fractional anisotropy at I = 3000 is shown in Fig- 
ures Q] and as a function of observed wavelength. The 
shaded regions are bounded by our two fiducial models 
(Population II stars with a Salpeter mass function will 
give the upper limit of the shaded region, while Pop- 
ulation III stars with a Larson mass function will give 
the lower limit). Other reasonable models, varying the 
mass or metallicity of the stars, will lie between these two 
limiting cases, since the amplitude of t he angular power 
spect rum will lie between these cases (jFernandez et al.l 
I2010D . We also show a range of f esc , from f esc = 1, 
(where all the ionizing photons escape from the halo into 
the IGM), to f esc = 0.1. Results are shown for a case 
where reionization progresses with a high efficiency, the 
minimum halo mass is 10 8 A/q, and small halos can be 
suppressed. However, these assumptions do not greatly 
affect the results. 

As seen in this figure, the escape fraction has a large 
effect on SI/I. The mean level of the CIB, given in equa- 
tion 21 has no dependence on the escape fraction. The 
angular power spectrum, given by equation [71 has a fac- 
tor of (1 — fesc)- Therefore, when f esc rises, the level of 
the nebular contribution to the angular power spectrum 
will fall. This causes the overall level of SI /I to fall. 

This drop-off of SI jl for larger values of the escape 
fraction is more pronounced at longer wavelengths. To 
see why this occurs, we can look at the mean spectrum of 
starlight and nebular emission for a high redshift galaxy 
in Figure [3] (The definition of the bands shown are given 
in Table [3J) In the near-infrared bands (A < 4 fj,m), 
there is always a large contribution from the stellar black- 
body emission. At longer wavelengths, the stellar emis- 
sion drops off very quickly, while, if the escape fraction 
is low, the nebular emission remains relatively high. If 
the escape fraction is high, however, the nebular emis- 
sion component of the angular power spectrum would 
be diminished. This is particularly noticeable at longer 
wavelengths. 

For small values of f esc , the range of allowed values 
for SI 1 1 is quite narrow. This range widens if f esc = 1. 
This is a consequence of the spectral shape for our two 
stellar models. For Population III stars, the blackbody 
spectrum of the star is steeper and there are more ioniz- 
ing photons to be processed into nebular emission. The 
stellar spectrum is almost always equal to or less than 
that of the nebular component. On the other hand, the 
stellar blackbody is greater than the nebular component 
at short wavelengths for metal poor Population II stars. 
As the wavelength increases, the emission from the stel- 
lar component will drop below the level of emission from 



Band 


llon+fal nf Rann [ /(in l 


V V OiVtiUdllVA 1 LL lib] 


AA/A 


T 

j 


1 25 


11-14 


9/1% 








1S% 


K 


2.2 


2.0 - 2.4 


18% 

J-O /o 


L 


3.5 


3.0 - 4.0 


99% 

/o 


M 


4.8 


4.6 - 5.0 


8.3% 


N 


11 


7.5 - 14.5 


64% 


Q 


21 


17 - 25 


38% 


z 


34 


28 - 40 


35% 


HI 


75 


60 - 90 


40% 


H2 


110 


90 - 130 


36% 


H3 


170 


130 - 210 


47% 



TABLE 3 

Band definitions used for infrared bands. HI, H2, and H3 

DENOTE THE BANDS OF THE HERSCHEL PACS INSTRUMENT. ALL 
BANDS ARE ASSUMED TO BE RECTANGULAR. 



the nebula. If f esc is small, the total emission of the ha- 
los for Population III and Population II stars is similar, 
so the range of SI /I is narrow. If f esc is large, the ampli- 
tude of the angular power spectrum of large Population 
III stars will be more affected than smaller Population II 
stars, widening the range of allowed values of SI /I. 

In Figure [2j we see the the fractional anisotropy for the 
case when galaxies contain dust. Here, we no longer see 
a decrease in SI /I at long wavelengths. This is because 
dust will reprocess the stellar and nebular emission, re- 
emitting the light at longer wavelengths. (An illustration 
of the SED with dust included is also shown in Figure [3] ) 
Unlike nebular emission, the dust component will not fall 
to zero when f esc = 1. In fact, the change between the 
dust emission from f esc = and f esc = 1 is only slight. 
Therefore, while the nebular component will disappear 
if fesc = 1) the dust component will still be present, 
causing the angular power spectrum at long wavelengths 
to remain high, raising SI /I. This is a direct result from 
the fact that nebular emission only results from ionizing 
photons, while photons of lower energies can be converted 
into dust emission. 

6. COMPARISON TO OBSERVATIONS 

Measurements of the CIB are notoriously hard to per- 
form. Finding an accurate mean is particularly difficult 
because precise foreground subtraction is needed. How- 
ever, observations, such as with CIBER, AKARI, and 
Herschel, continue to improve, resolving foregrounds in 
more detail and obtaining more reliable observations for 
the CIB. 

Many observations have been made to try to under- 
stand the contribution of high redshifts to the CIB from 
1-4 microns. In order to uncover any residual emission in 
the mean or fluctuation observations, one must carefully 
take into account all of the foreground components. 
Zodiacal light is a major contaminant, and because it 
is very difficult to model, it is not straightforward to 
subtract from the CIB. In addition, foreground galaxies 
at lower redshifts must be taken into account. Despite 
the difficulty, there have been many attempts to measure 
the mean level of the CIB in the near-infrared due to 



high r edshift stars dDwek &; Are ndt 199 81 iGoriian et al 



200C 


; iKashlinskv & Odenwald 120001: IWrieht & Reese 


200C 


: Wright 20011: iCambresv et al. 


20011: Totani et al. 


2001 
Kasl 


: IKashlinskv et al.l 
llinskv et al.l 12007a 


2002L2004 
IhL 120121 


: IKashlinskv! 2005: 
Maeliocchetti et al. 



THE COLOR OF FRACTIONAL ANISOTROPY 



■5 




1 10 100 

Wavelength [microns] 



Fig. 1. — The Fractional Anisotropy, SI/I, as a function of wavelength, for various values of f eS c, in cases without dust. Population II 
stars with a Salpeter mass spectrum provide the upper limit of the shaded regions, while Population III stars with a Larson mass spectrum 
provide the lower limits. Other reasonable assumptions for the mass and mctallicity of the stellar populations should lie within the shaded 
regions. 




10 

Wavelength [microns] 



100 



Fig. 2. — The Fractional Anisotropy, SI/I, as a function of wavelength, for cases with dust. 



200l lOdenwald et all l200l ICoorav et al l [200l 
Matsumoto et al J l2005t iThompson et alJl2007al lb1). 

Fluctuation observations are, in theory, easier to per- 
form, since they do not need an accurate zero point, 
and instead rely on variations from one region of the 
sky to another. However, these observations still rely on 
careful and complete subtractio n of foreground sources , 



and a l so remain controversial ( Kashlinskv et al 



20?)7bl l20H ICoorav et al1 l2007l: 



Thompson et al.l 



[2005 



Matsumoto et alJl2011j )T 

Observations are even more difficult in the mid 



2007at 



and 



far-infrared. One problem is that foregrounds that were 
present in the near-infrared are even more prevalent in 
the mid and far-infrarcd. Zodiacal light peaks at about 
20 \xm, which washes out most detections of the CIB in 
this range. In addition, Galactic cirrus is a main contam- 
inant, however, observations arc possible in clean regions 
of the sky. Finally, as wavelength increases, foreground 
galaxies become more difficult to resolve. All of these 
problems lead to only a fraction of the CIB in the mid 
and far-infrared being resolved into low-redshift galax- 
ies. It is likely that only a very small (and currently 
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Fig. 3. — The spectra of a 10 9 Mq galaxy at z 10. We assume a star formation efficiency of /* = 0.1. Two stellar populations are 
shown - Population II stars with a Salpeter mass spectrum, or Population III stars with a Larson mass spectrum. In addition, two limiting 
cases of the escape fraction are shown - f eac = and f esc = 1. When f esc = 0, the ionizing radiation stays within the galaxy, creating 
nebular emission, which is not present in the halo if f esc = 1. The dust emission will fall only slightly when f esc = 1. Bands are denoted 
by the shaded gray regions. Hershel PACS bands are labeled as HI, H2, and H3. 



unknown) percentage of this excess is from z > 6, so 
care must be taken in interpreting observations. 

There has been a great push to understand the 
CIB at longer w avelengths. At 100 and 160 fj,m, 
iPenin et al.l ( 20111 ) measured the mean and fluctuation 
power of galaxies at all redshifts using observations 
from IRIS/IRAS and Spitzer/MIPS, respectively. 
Fluctuations of the cumulative CIB have been taken 
in the mid-infrared to submi llim eter wavelengths 
(iKashlinskv k Odenwaldl [20001 lLagache k Pueetl 
2000J iMiville-Deschenes et al.l 120021 iGrossan k Smood 
20071 : lLagache et all 120071; lAmblard et al.l l20lTF 



Matsuura et all 120111: 
20111: IPvo et"all 12012ft . 



Planck Collaboration et all 
The cumulative mean level 
of the CIB from g alaxies at al l redsh i fts has been 
measured as well dFixsen et al 119981 Mauser etiil 
19981: lLagache et all 120001: iWrightl 12001 lOdegard et al 



2007 



Matsuura et all 120111) . The mea n CIB has 



been measured as a function of redshift dBerta et all 
120111: Uauzac et all 120111: iBethermin et all 12012ft . while 
measurements from BLAST and Planck from 250 
to 1400 microns (jPlanck Collaboration et all 12011ft 
could indicate that galaxies at a higher redshift (here, 
z > 1.2 — 2) could contribute more to the CIB as the 
wavelength increases. Curre n tly, th e best measur e ments 
(jPlanck Collaboration et"all I20TTI: iViero et all 12009ft 
show that the fractional anisotropy is at the order of 
15%, however, these measurements include galaxies at 
all redshifts. 

One way to subtract unresolved low redshift galaxies 
in the mid to far-infrared to a more complete level is 
to use a stacking algorithm. This typically involves us- 
ing the locations of known galaxies at a shorter wave- 
length, stacking these locations of a longer wavelength 
image, and utilizing this stack to calculate the CIB ac- 



counted from these galaxies. If stacking is relied upon, 
more of the CIB at long wavelengths c an be resolved 
into lo wer redshift galaxies. For example. iMarsden et all 
(|2009ft used stackin g to resolve 100% o f the CIB as de- 
tected with FIRAS (iFixsen et aflll998ft at 250 350, and 
500 microns using BLAST. iDole et all (|2006l ) used 24 
jum sources from Spitzer/MIPS data to stack images at 
70 and 160 /im. They were able to resolve 79%, 92 %, 
and 69% of the C IB at 24, 70, and 160 \im respectively 
iBerta et all (|2010j ) resolved 45% and 52% (without stack- 
ing) and 50% and 75% (with stacking) of the CIB at 100 
and 160 am using He r schel/ PACS data. At longer wave- 
lengths, iGreve et all (f20ToTi resolved 16.5% of the CIB 
at 870 microns using stacking. While it is possible that 
some of the remaining flux is from low redshift galax- 
ies, it is also possible that some of this unresolved CIB 
could be due to high re dshift galaxies. (See, for example, 
iMatsuura etall (|2011ft ). 

We compare some of these observations to our mod- 
els for the fractional anisotropy at high redshifts. While 
precise measurements of the mean are challenging, mea- 
surements of the fluctuation power are becoming more 
reliable. Because of this, in Figure [4] we show the frac- 
tional anisotropy predicted using various recent observa- 
tions of the fluctuation power, assuming an upper limit 
of the mean CIB due to high redshift stars is either 10 
nW m~ 2 sr _1 or 1 nW m~ 2 sr _1 . Because it is unlikely 
that a z > 6 component of the mean CIB will be much 
higher than this, very low values of SI /I, and thus very 
high values of the escape fraction from a dust-free popu- 
lation, can be ruled out. More definitive conclusions can 
be reached as observations continue to improve and our 
understanding of foreground emission grows. 

It is important to remember that there is still 
low-redshift contamination contributing to both the 
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Fig. 4. — The Fractional Anisotropy, SI/I, in comparison to recent observations. The upper limit of the mean CIB from z > 6 stars is 
assumed to be either 10 nW m — 2 sr -1 (lower set of arrows) or 1 nW m -2 sr _1 (upper set of arrows). These assumptions provide the lower 
limits of SI/I. Shown are observations from Kashlinsky et al (2012) (red arrows), Matsumoto et al (2011) (purple arrows), Cooray et al 
(2007) (green arrows), Pyo et al (2012) (blue arrows), and Matsuura et al (2011) (pink arrows). The horizontal hatched region shows cases 
with dust and f eBC = 1, the diagonal hatched region shows cases with no dust and f eBC = 1, the grey shaded region shows cases with dust 
and fesc = 0.1, and the black region shows cases with no dust and f eS c = 0.1. 



mean and fluctuations, especially at longer wavelengths. 
Therefore, these results should be interpreted with this 
in mind. As observations improve, these results will be- 
come more reliable. 

7. CONCLUSIONS 

We have shown the observable signatures of high 
redshift populations with different values of the escape 
fraction f esc and dust content in the CIB. It is possible 
to distinguish these populations through observations of 
the fractional anisotropy of the CIB. The global escape 
fraction of high redshift galaxies is a main variable that 
can be probed in this way, since the angular power 
spectrum is dependent on it, while the mean is not. In 
addition, dust will transform the SED of the galaxy, 
thus leaving an imprint on the fractional anisotropy. 
Therefore, low values of the fractional anisotropy will 
be indicative of a population of stars with a high escape 
fraction and little dust. This will be more noticeable at 



longer wavelengths. While observations are still difficult, 
improved observations could be able to distinguish 
between these populations. 
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